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of conservation of energy balance is presented, depends on its type of heat transfer within a whole 
system of soil and foundation. In order to understand the mechanical implication behind their application, 
simulation numeric with finite difference method is conducted, concerning an energy pile in homogenous 
soil under static thermal load. The study takes into account two different conditions of contact between 


sch dati soil and pile: perfectly contact and sliding contact using frictional interface elements. The results show 
Building energy that temperature-induced mechanical behaviour of pile and soil is strongly related to the condition of 
Heat transfer contact between them. Further work is projected to consider a more appropriate law that corresponds to 
Thermal load cyclic thermal loading of energy piles due to its seasonal cooling and heating operation throughout the 
Thermo-mechanic behaviour year. 

Condition of contact © 2012 Elsevier Ltd. All rights reserved. 
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the average ground temperature at shallow depth 10-50 m, where 
its value ranges from 10°C to 15°C in most European countries 
[1]. This kind of technology is well known as conventional ground 
heat exchanger, where ground thermal energy is extracted from the 
soil via heat pump connection to the building. The ground operates 
as heat source supplying warm energy to building during winter 
season, whilst during summer season it functions as heat sink when 


E-mail addresses: msuryatr@polytech-lille.fr (M.E. Suryatriyastuti), a o of building is abe P l dh h 7 
hussein.mroueh@polytech-lille.fr (H. Mroueh), sebastien.burlon@ifsttar.fr There are two types of conventional ground heat exchangers. 


(S. Burlon). The first one is open loop system where the groundwater is used 


1. Introduction 


Utilisation of shallow geothermal energy in providing thermal 
needs of building is already common. Shallow geothermal energy is 
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Fig. 1. Energy pile components. 


as a heat carrier and pumped directly from an aquifer to the heat 
pump. The remained waste water is returned to the same aquifer 
or to another well. This system is simpler but it induces some envi- 
ronment problems like groundwater pollution and needs high cost 
of maintenance due to clogging in the wells and heat exchangers 
components. 

The second one is closed loop system where closed-coils 
absorber pipes are embedded into borehole, laid either horizontally 
or vertically. A heat carrier medium is circulated through absorber 
pipes to transport heat from the ground to the heat pump. This sys- 
tem is also called as borehole heat exchanger (BHE) system. Because 
of its high cost of drilling (for vertical loop) and its need of huge land 
area (for horizontal loop), closed loop system has recently been 
merged into structural foundation elements. 

Energy foundations have both functions as structural elements 
and heat exchanger elements. This system is like closed-loop con- 
ventional ground heat exchanger system unless the absorber pipes 
are installed on foundation elements. The main advantage of this 
system is the foundation elements that serve as heat exchang- 
ers are already required for structural reason and no need to 
be constructed separately. Furthermore, concrete as heat transfer 
medium in energy foundations produces higher thermal energy 
due to its higher thermal capacity than grout or rock in bore- 
hole heat exchanger. Therefore, as well as reduce significantly cost 
installation and land use, energy foundations are able to increase 
thermal productivity while preserving environmental sustainabil- 
ity. 

Some applications of energy foundations technology have been 
developed in last decade, such as: Main Tower in Frankfurt- 
Germany on 1999(112 energy piles) [2], Keble College in Oxford-UK 
on 2001 (90 energy piles, 74 MWh annual heating and 55 MWh 
annual cooling) [3], Dock Midfield Zurich Terminal Airport- 
Switzerland on 2003 (306 energy piles, 2210 MWh annual heating 
and 1170 MWh annual cooling) [4], and Lainzer Tunnel Vienna- 
Austria on 2004 (59 energy piles, 214 MWh annual heating) [1]. All 
of those projects confirm that utilisation of energy foundation sys- 
tem can save the annual overall costs, even though in the beginning 
the investment costs are highly expensive. Energy piles applica- 
tion at Zurich Airport reported cost saving up to €36,900 relative 
to conventional system with annual thermal cost €0.04/kWh [4]. 

Despite of its successful implementation, there are two main 
problems of this technology to be considered. One concerns about 


sustainability of thermal capacity in the system, another concerns 
about mechanical durability of structure foundation due to thermal 
solicitation [5]. Derived from those main problems, this paper will 
discuss separately each problem of particularly energy piles into: 


(a) physical process of heat transfer which is influenced by seasonal 
ground temperature and groundwater flows (Section 3), 

(b) thermal transfer balance between soil-concrete foundations- 
and building (Section 4), 

(c) mechanical implication induced by additional thermal loads 
which change the behaviour of structure (Section 5). 


In order to understand energy piles’ thermo-mechanical 
behaviour, preliminary study of an energy pile under static ther- 
mal loads in each operation season is being simulated. Since the 
response of piles to vertical loads strongly depends on the char- 
acteristics at the soil—pile interface, it is of major interest to take 
into consideration the interface condition in the analysis of the 
behaviour of piles under thermal loading [6,7]. In this study, two 
comparisons of contact conditions are proposed: perfect contact 
and sliding contact between soil and pile using frictional interface 
elements (Section 6). 


2. Energy piles system 


Energy piles system contains of three principal circuits (Fig. 1). 
The first circuit is the soil as heat source in winter season where 
ground energy is withdrawn for heating needs of building. It 
becomes inversely heat sink in summer season where heat is 
rejected back to soil in order to recharge ground energy and to cool 
the building. 

The second circuit contains of closed absorber pipes in concrete 
piles through which heat carrier fluids are circulated to exchange 
energy from/to soil to/from building. Heat carrier fluids are water 
with antifreeze solution which is able to avoid freezing at incli- 
nation angle of absorber pipes. Hence, glycol-water mixtures are 
the most suitable solution, containing also additives to prevent 
corrosion in the header block of valves of the heat pump [8]. Con- 
cerning the type of absorber pipes used in energy piles, Hamada 
et al. [9] conducted performance test of three kinds closed absorber 
pipes: U-shaped, double U-shaped, and indirect double-pipe in the 
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building project constructed at Sapporo, Japan. As a result, the U- 
shaped pipe type was employed from the viewpoint of economic 
efficiency and workability due to its smallest heat rejection rate [9]. 

The third circuit is closed-pipelines heating/cooling network 
which embedded in the floors and walls of building. Heat pumps 
are linked between second and third circuit in order to raise the 
original non-usable heat resources to a higher usable room temper- 
ature. Quantity of energy production depends on the efficiency of 
heat pumps, which is described by the coefficient of performance 
(COP). COP is the ratio of the output energy divided by the input 
energy (electricity for the compressor) and varies from 3.0 to 6.0 
with present equipment (the higher the number the better the effi- 
ciency) [10]. For a best energy produced with economical cost, heat 
pumps require a number of COP higher than 4.0 [8]. 

Two operation schemes of energy piles are possible: single mode 
of energy extraction/energy input and seasonal mode operation 
with heating and cooling storage (Fig. 2). For single mode opera- 
tion, the energy flow takes place in one direction only, for example: 
heating purposes in glacial countries or cooling purposes in tropical 
countries. If only heating or cooling is performed, high perme- 
ability ground and high hydraulic gradient of groundwater are of 
advantage [8]. However, the seasonal operation uses the thermo- 
dynamic inertia of the soil in order to store thermal energy in 
the ground for later operation with reversed energy flow. Conse- 
quently, the seasonal operation can produce energy equilibrium 
automatically in the ground over acomplete heating/cooling period 
of a year [1]. This operation system is preferred by the ecologist 
due to its environmental friendly for the sustainability of ground- 
water. 
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Fig. 2. Scheme of energy piles system operation: (a) energy extraction and (b) energy storage. 


Heat transfer in energy piles system is surrounded by complex 
process. In the first circuit, the transfer depends on soil ther- 
mal and hydrology properties, which are specifically conducted by 
thermal conductivity, specific heat extraction, and water content 
of soil. In the second circuit, thermal transfer affects mechanical 
behaviour at concrete piles, predominantly related to the inter- 
action between soil-concrete piles-and upper structure. Finally, 
coefficient of performance of heat pumps is the key to avoid 
the energy loss between the second and third circuit. Because of 
the fact that no design method is yet available to consider the 
complex interaction between thermal storage and the mechanical 
behaviour of energy piles [11], the present work offers some steps 
to implement energy piles’ technology as explained at flowchart in 
Fig. 3. 


3. Physical process in heat transfer 


Soil is a porous medium formed by solid grains, water, and 
gases. Due to their compositions, heat transfer in soil consist of 
many mechanism, including: conduction through the soil grains, 
liquid, and gases; convection through liquid and vapour diffusion; 
latent heat transfer through evapouration-condensation cycles 
during the phase change of water; radiation in gas-filled pores; 
and freezing-thawing process that should be avoided in energy 
foundations system [8,12]. Heat conduction leads dominantly, heat 
convection in liquid diffusion occurs by presence of groundwater 
flow, vapour diffusion and latent heat transfer might only possi- 
ble in unsaturated soil, and radiation contributes negligibly to heat 
transfer. 
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3.1. Ground temperature 


Burger et al. [13] divides the soil into three zones of temperature: 
heterothermal, neutral, and homothermal zone. External temper- 
ature and radiation affect ground temperature in heterothermal 
zone. In this zone, temperature varies in time and depth as function 
T(z,t): where T(t) defines thermal cycle loads in annual exter- 
nal temperature and T(z) indicates thermal waves’ propagation in 
depth. This is the most critical zone depth that should be properly 
taken in account of energy piles design calculation. The function of 
temperature forms a sinusoidal wave, as following equation: 
T(z, t) = Tave + Ape -7/9 (sin (ot = =)) (1) 
where Tave is the average annual soil temperature; Ap is maximum 
annual amplitude; z is the depth in which temperature is investi- 
gated; d is the damping depth of annual fluctuation; and w is the 
annual radial frequency. 

Fluctuation of temperature becomes constant since it arrives 
in neutral zone at depth 10m until 50m below the surface [13] 
and remains closely at average air temperature, where the value 
depends on altitude h: 


Tair = 11.3 — 0.52h (2) 


In homothermal zone, temperature is no longer influenced by 
external thermal flux and no longer varies by annual sinusoidal 
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Fig. 3. Flow chart of design steps. 


waves. This zone is where the deep geothermal or internal terres- 
trial flux takes place. Temperature increases regularly in depth with 
geothermal gradient 3 °C per 100m [14]. 


3.2. Soil thermal properties 


Heat flux propagation in soil is transmitted by heat capacity and 
thermal conductivity. Volumetric heat capacity C, is a product of 
mass volume p with specific heat extraction c that indicates the 
ability of a substance to store heat energy. The greater its heat 
capacity, the more heat it can gain or lose per unit rise or fall in 
temperature. Whereas thermal conductivity à is the quantity of 
heat transferred through a unit area of the conducting body in unit 
time under a unit temperature gradient [15]. This indicates the rate 
at which heat will be transferred between exchanger loops and 
the surrounding transfer medium for a given temperature gradi- 
ent. Thermal conductivity of the surrounding transfer medium is 
the critical value that determines the number of heat exchanger 
pipes required to produce energy needs of building by pumping 
specific fluid temperature through the ground loops. 

Thermal properties of soil, including heat capacity and thermal 
conductivity, cannot be separated with mineral composition of soil 
itself as mentioned in Eqs. (3) and (4) with index “s” denotes solid 
composition and “f” denotes fluid composition. Porosity factor n, 
water and air content will lead the quantity of heat produced. Moist 
or saturated soils have greater thermal capacity than dry soils, since 
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thermal conductivity for dry soil is about five times lower than that 
of a saturated soil [14]. 


Cy = (1 — n)pscs + npr Cp (3) 
A=(1—n)As + nay (4) 


3.3. Groundwater flow 


As water and air content take important place in producing ther- 
mal capacity of soil, the presence of groundwater flow is one of key 
factor in heat transfer mechanism. In fully saturated soil with no 
presence of pore gas, groundwater composes only of liquid flow, 
so that heat convection occurs in liquid diffusion. In unsaturated 
soil where the pore gas presences inevitable, groundwater flow has 
both of liquid and vapour flow, so that heat convection counts both 
liquid and vapour diffusion. It is also possible to have a water phase 
change into fully vapour that correlates with vapour latent transfer. 
A widely used theory for coupled heat and moisture flow through 
soils was developed by Philip and de Vries [16] for both liquid and 
vapour phase flows. 

Groundwater flow in soil varies under water content 6 and tem- 
perature T gradients. Water content variation is led by isothermal 
diffusivity Dg while temperature variation is controlled by thermal 
diffusivity Dr. In general, water content gradient is less than those 
of temperature, therefore the flux due to water content gradient 
is more important than those due to temperature gradient, even 
with the maximum temperature gradient [16]. Hadley and Eisen- 
stadt [17] have observed that moisture transfer under temperature 
gradient is negligibly small both in very dry (unsaturated) and in 
very wet (saturated) media [16,17]. 

In liquid phase, the velocity of groundwater is given by Eqs. 
(5)-(7), referring to Darcy’s equation where k is hydraulic con- 
ductivity and iis unit vector in vertical direction. Here, isothermal 
liquid diffusivity Dg; and thermal liquid diffusivity Dr, depend on 
capillary suction head w and surface capillary tension o. 


Ù, =—De grad 0 -— Dr grad T- ki (5) 
_ av 

Do. = kag (6) 
_ Vdo 

Drı = kg aT (7) 


Movement in vapour phase is a process of the diffusion of water 
vapour in the air-filled pores which can be approximated by mod- 
ifying the Fick’s law of diffusion: 


Ùy = -Dr „grad T — Do „grad 0 (8) 


where Dr, is the thermal vapour diffusivity and Dg , is the isother- 
mal vapour diffusivity. 

Groundwater flow is then a sum of flow in liquid and vapour 
phase. The two liquid diffusivities tend to be the most important 
ones at high moisture contents, whilst the two vapour diffusivities 
are dominant at low moisture contents [16]. Refers to precedent 
hypotheses, groundwater flow in saturated soil is given in Eq. (9) 
while in unsaturated soil is mentioned in Eq. (10). 


Vsat = 4 grad 0=-—k grad w (9) 


Vunsat = —Dg „grad 0 (10) 


Site installation of Finkernweg at Switzerland has been con- 
ducted by Fromentin et al. [5], concerning 75 energy piles with 11 m 
active length to produce 90 kW and 7000 m? heat energy. Thermal 
measurement is investigated during February until May 1995. It is 
shown that groundwater flow does not have significant influence 
in short term operation. However, groundwater flow may repro- 
duce heat concentration energy around the piles with velocity up 


to 10cm/day in long-term period [5] so that it plays an important 
role for pile’s dimensioning. Also, if a flow of water is present and 
sufficiently large, a natural regeneration of the soil is achieved [14]. 
As aresult, the amount of heat extracted during winter (for heating) 
no longer depends on the amount injected in the summer (for the 
cooling). Hence, groundwater flow is of advantage to single mode 
energy extraction system but further study is required to ensure 
sustainability of ground water. 


3.4. Heat flux 


The rate of energy transfer in soil through solid, liquid, or gas 
pores is called as heat flux. It is important to define the heat flux in 
order to quantify the energy produced. 

For heat conduction mechanism, heat flux is based on Fourier 
law: 


cond =-A grad T (11) 
while in liquid convection mechanism is determined by: 
conv = prcfUi(T — To) (12) 


It should be noted that groundwater flow in heat convection 
flux is the groundwater flow in liquid diffusion. In reverse when 
we consider latent heat flux, groundwater flow which contributes 
in energy transfer is groundwater flow in vapour diffusion, as men- 
tioned as: 


latent = Lydy (13) 
3.5. Energy conservation 


Finally, heat transfer in soil should satisfy the energy conserva- 
tion of thermodynamics law, which is: 


oT 
ot 
Table 1 summarises heat transfer mechanism for each type of 


soils and their energy conservations are given in Eqs. (15)-(18). 
For saturated soil with no groundwater flow: 


G=+Vqg=0 (14) 


a =div(rA grad T), (15) 
for saturated soil with groundwater flow: 
ea —prcrk grad y grad T=div(A grad T), (16) 


for unsaturated soil without groundwater flow: 


Coan = dinà grad T)—Ly div(Do vap grad 6), (17) 


and for unsaturated soil with groundwater flow: 


Coan — perk grad w grad T = div(Àà grad T) 


-Ly div(Devap grad 0) (18) 


4. Thermal transfer balance 


Energy piles are dual-purpose as structural foundation and heat 
exchanger. Thermal capacity of their heat exchanger elements 
will take a role for structural dimensioning of their foundation 
elements. Their functions to satisfy thermal and mechanical (struc- 
tural) capacity should be carefully designed by mechanical and 
geotechnical engineers. For Q3 is the energy required of building at 
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Table 1 
Division of heat transfer process in soils. 


Heat transfer Saturated 


Unsaturated 


No ground-water flow 


With ground-water flow 


No ground-water flow With ground-water flow 


In solid conduction v v 
In liquid convection x 
In vapour latent x x 


< 


v 
X v 
v v 


< 


Table 2 
Properties of model. 


Thermo-elastic properties 


Unit Soil Concrete 


Density 

Bulk modulus 

Shear modulus 

Thermal conductivity 

Specific heat extraction 
Coefficient of thermal expansion 


RO Framed 


Thermo-elastoplastic properties Unit 


Pa 1950 2500 
MPa 10 
MPa 3.75 7500 
W/m? 1.5 1.8 
J/kg°C 800 880 
ime 5 x 10-6 1x 10-5 


Interface 


Normal stiffness 
Shear stiffness 
Cohesion 
Friction angle 


~ 
3 


MN/m 
MN/m 
kPa 


~ 
s 


ea 


22 
8.33 
1 

30 


third circuit and Qə is input energy at heat pump at second circuit, 
the energy conservation which occurs in the heat pump is: 


Q3 = Qo + Pup (19) 
; 1 : 
@ (1- ap) = % (20) 


Input energy at heat pump Qz is equal to resources heat energy 
from soil through concrete piles Q4. Hence, to design the number 
n of concrete piles required, Eqs. (21) and (22) are given as below, 
where | is the length of material. 


Qo = Qı (21) 
J Arddio = [CC (22) 


To minimise energy loss between each circuit, high-raised 
energy can be produced by the help of heat pumps performance. 
Study of thermal performance of energy piles applied at Zurich Air- 
port showed that annual heat pump performance coefficient is 3.9. 
The overall system efficiency, defined by the ratio of the thermal 
energy delivered by the system (heating and cooling) over the total 
electric energy required to run it, is established to 5.1 [4]. With 
a high efficiency system, energy piles will maintain in long-term 
operation with less energy loss. 


5. Mechanical equilibrium 


In general, the mechanical response of a system must sat- 
isfy equilibrium, compatibility, material constitutive behavior and 
boundaries conditions [18]. Among these conditions, thermal effect 
is involved in mechanical response of energy piles. Indeed, the 
change of temperature both in concrete and soil will produce 
supplementary thermal deformation that will modify the strain 
state of materials. The quantity of thermal deformation eth is a 
product of thermal expansion coefficient œ with gradient temper- 
ature occurred as mentioned in Eq. (23). Elastic deformation ¢° 
will change in function of thermal deformation obtained and total 
deformation € caused by static weights of system structure, given 
in Eq. (24). 


eth = «AT (23) 


e =e-eth (24) 


Since lateral movements of concrete pile are restrained by soils, 
thermal deformation will mainly occurs in vertical axis of pile. The 
vertical displacements are not uniform in depth and strongly relates 
to lateral friction mobilised in the shaft of pile, in contact with soil 
[19]. The supplementary thermal deformation brings out increas- 
ing of mobilised shaft friction at interface soil-concrete piles also 
uplift displacement and additional thermal compressive stress at 
interface concrete piles—upper structure [20]. The relation between 
axial force Q, mobilised shaft friction qs, and deformation at a given 
depth z is expressed on following equation [21]: 


Q(z) = €°(Z)AEpite (25) 


Q(z) _ AEpite Ae*(z) 
9s(2) = 4 zy = A(z) dz 


(26) 


where Epil is young modulus of concrete pile, A is area of pile’s 
section, As is shaft area at certain depth. 

A four storey building built with 97 piles was carried out 
to an in situ test of a heat exchanger pile at the Swiss Federal 
Institute of Technology in Lausanne (EPFL). Data from tests indi- 
cate that a temperature increment of 1°C results in an additional 
temperature-induced vertical force on the order of 100 kN [22]. As 
a consequence, the total axial load in the pile is twice as large as the 
one due to purely mechanical loading, with a large solicitation of 
the toe because the pile toe sits on very stiff sandstone that restrain 
the pile. 

Another in situ test of energy piles was conducted in Lambeth 
College, United Kingdom with a total 143 piles. The different with 
those from Lausanne test is in this test, the piles are floating in the 
ground so that both ends of pile are free to move. When a cooling 
cycle is applied to the pile, the pile contracts and leads to induce 
tensile force. Along the upper part of the shaft, mobilised lateral 
friction on the interface soil-pile will increase in the same sense 
as that mobilised by compression loading applied at the pile head, 
while in the lower part of the pile will be in the opposite sense 
so reduce the mobilised friction [23]. The reverse responses are 
appeared when a heating cycle is applied, which suggest the pile 
behaves in purely thermo-elastic. 

If the energy piles work on seasonal mode operation sys- 
tem, they are subjected by alternate thermal loads. The alternate 
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Fig. 4. 3D mesh discretisation of model: (a) one-fourth of model and (b) interfaces location at Model 2. 


contraction/expansion piles create cyclic behaviour at interface 
soil-concrete piles-and upper structure. As explained in Section 
2, from an energy balance point of view, seasonal mode operation 
system is the most favourable one referring to energy balance point 
of view. But if we look from the view of mechanical durability, 
cyclic behaviour at interface will decrease strength resistance of 
structure. 


6. Case study: temperature-induced mechanical behaviour 
of an energy pile under static thermal loading 


6.1. Numerical model 


Numerical method such as finite difference method provides an 
efficient tool for solving complicated partial differential equations 
such as those that govern the thermo-hydro-mechanical problem 
of energy pile, because of its ability to integrate complexities aris- 
ing from: geometry, thermo-hydro-mechanic loading, interaction 
between different rigidity bodies, etc. The case under study con- 
cerns with the analysis of an energy pile under static thermal 
loading. A single pile with a square section of B = 60 cm and length 
L=15m is founded into a homogeneous soil mass. Both of pile and 
soil are assumed to behave into linear thermo-elastic domain. Prop- 
erties of materials are shown in Table 2, where thermal properties 
are chosen in reference to Hillel [15]. 

Analyses are carried out using finite difference code FLAC3P 
[24]. Fig. 4a shows the numerical mesh used into the numerical 
modelling. Owing to the symmetry of the problem (symmetry of 
load and geometry), only one-fourth of the complete domain is 
modelled. The mesh was fixed after a parametric analysis con- 
cerning the lateral extension of the soil domain. The soil is refined 
around the pile, in order to increase precision in the areas of high 
strain gradient. According to Karthigeyan et al. [25] and Mroueh 
and Shahrour [7], the soil mass dimension depends on the pile 


section and length. Lateral extension of the domain is equal to 15 m 
(25B). The height of soil mass is equal to 30 m (2 L). As mentioned 
in Section 1, the response of piles to vertical loads strongly depends 
on the characteristics at the soil-pile interface [6,7]. This study 
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Fig. 5. Fluctuation of initial ground temperature. 
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Fig. 6. Initial condition of ground temperature. 


proposes comparison of two conditions of soil—pile contact in 
order to understand their different behaviour. Model 1 considers 
a perfect contact between soil and energy pile that means there is 
no relative movement between the nodes of pile and soil. Model 2 
introduces frictional model for the soil—-pile interface, governed by 
an elastic perfectly plastic behaviour. These interfaces elements 
are able to localise relative movements and frictions mobilised 
between two different solids (Fig. 4b). Elasticity part is described 
by the following relation: 


a J [s e (i) (27) 


where T and on denote the shear and normal stress and £t and €En 
indicate the shear and normal strain at the interface, respectively. 
The interface shear stiffness ks and the normal stiffness kn are 
chosen in accordance with interfaces’ theory and background of 
the code FLAC?) [24]. Plasticity part is defined by a yield surface 
based on the Coulomb failure criterion: 


[tT] =0n tan gj +c (28) 


where g; and c; are the friction angle of interface and cohesion 
factor at interface. The interfaces properties are shown in Table 2. 

At first, soil system is initialised by sinusoidal surface tempera- 
ture with average ground temperature 14°C. The soil has internal 
terrestrial heat volume r as 0.001 W/m? and terrestrial surface flux 
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Fig. 7. Normal forces due to gravity load. 


qat depth 30 mas 0.0544 W/m?. Due to symmetrical condition, flux 
null is imposed in symmetrical axes. The presence of groundwater 
flow is neglected in assuming that soil is fully saturated so that heat 
diffusion works in purely conduction. This sinusoidal surface tem- 
perature produces fluctuation of ground temperature in depth and 
in time that becomes constant near to average ground temperature 
at the depth >10 m (Fig. 5). This is justified by theoretical hypothesis 
in Section 3.1. 

According to that fluctuation, seasonal inlet fluid temperature is 
applied into the pile and assumed to be uniform along the pile’s axis. 
The model is conducted in two separated operation time: (a) winter 
mode operation with inlet temperature is set on 5 °C (cooling cycle 
- Model W) and(b)summer mode operation with inlet temperature 
25°C (heating cycle - Model S). Initial temperature of soil before 
the injection of inlet fluid at each season is shown in Fig. 6. 

The pile is only subjected by gravity load and then by temper- 
ature applied during six months loading for each operation mode. 
No service load coming from upper structure is considered to prop- 
erly understand the effect of thermal load in pile’s mechanical 
behaviour. Fig. 7 shows the initial normal forces due to gravity load 
in two different conditions of contact. For the clarity, positive sign 
signifies compression on forces and settlement on displacements. 
Model 2 with interfaces imposed in the zone of contact gives reduc- 
tion of normal forces in pile because interface zones take a role 
to transfer the forces more uniformly to the system by reducing 
discontinuity of residual forces. 

Comparison of thermal diffusion by the function of contact’s 
condition shows a remarkable influence (Figs. 8 and 9). Within the 
presence of interface elements, there is high concentration of ther- 
mal field at surface contact between pile and soil. This leads to a 
statement that interface behaviour depends also on temperature as 
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(a) 


(b) 


Fig. 8. Field of temperature after cooling cycle: (a) Model 1 perfect contact and (b) Model 2 with interfaces. 


(a) 


(b) 


Fig. 9. Field of temperature after heating cycle: (a) Model 1 perfect contact and (b) Model 2 with interfaces. 


a function parameter, that interface nodes take a role on tempera- 
ture transfer. It is suggested to conduct further study to look about 
interface behaviour law with the coupling of temperature function. 

Concerning to the mechanical behaviour due to thermal loading, 
thermal displacement of the surface and along the pile are shown 
in Figs. 10 and 11. Winter mode operation produces a settlement 
at surface plus shortening of pile while summer mode operation 
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Fig. 10. Thermal surface displacements. 


leads to an uplift surface with expansion of pile. Contraction in 
winter mode induces tensile forces at the pile while dilatation in 
summer mode adds compressive forces at the pile. Model 2 with 
interfaces imposed gives clearly different mechanical responses 
and also reduces the value of temperature-induced normal forces 
produced (Fig. 12). 

Fig. 13 shows frictions mobilised along interface soil-pile at 
Model 2. Due to its contraction in winter mode, frictions mobilised 
along the upper part of pile increase whereas along the lower part 
of pile decrease. Inversely in summer mode, a reduction of frictions 
mobilised is appeared at upper part of pile while the lower part of 
pile has additional frictions mobilised. At middle part of pile, the 
friction mobilised tends to be constant as what it’s called null point. 
This behaviour is only occurred when both ends of pile are freely 
restrained. These results are in good agreement with in situ tests 
conducted at Lambeth College, UK [23]. 


6.2. Perspectives 


According to this preliminary numerical model, static thermal 
load in only cooling/heating operation give slight influence on 
mechanical behaviour of pile. Alternate contraction and dilatation 
during seasonal operation time is predicted to be a dominant factor 
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Fig. 11. Thermal vertical displacements of pile. 
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Fig. 12. Temperature-induced normal forces in the pile. 
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Fig. 13. Temperature-induced frictions mobilised at interface soil-pile. 


that provokes a strong degradation of structure’s resistance, espe- 
cially at interface soil-structure. 

Our perspective is to determine the interface behaviour between 
soil-structure due to thermal cyclic solicitation. To achieve this 
goal, it is better to consider an appropriate law that can define cyclic 
behaviour of interface soil-structure in taking into account plastic- 
ity condition. It should consider also the effect of delta temperature 
occurred at interface behaviour. On-going works have developed 
a constitutive law for interface soil-structure under mechanical 
cyclic loading, based on previous work of Shahrour and Rezaie [26]. 
This law takes into account kinematic hardening rules that allow 
to describe kinematic surface in each cycle loading. The further 
research aims to adopt this constitutive law into thermal regime, 
coupled with mechanical service loads. 


7. Conclusions 


Energy piles are bi-function foundations which support the 
loads of structure and provide thermal energy to the building 
as heat exchangers. They’re using shallow geothermal energy as 
heat source, successfully considered as an innovative environmen- 
tally friendly building’s technology. Energy piles can work in two 
system operations: (a) only cooling/heating operation system for 
glacial/tropical countries and (b) seasonal cooling-heating opera- 
tion system for four-season countries. The latter one is the most 
favourable for its automatic energy balance throughout the year. 

This paper covered a sum of literature study to identify the 
remaining problems in energy piles’ application. Quantity of heat 
energy produced depends on type of soil, saturation degree, pres- 
ence of groundwater flow, and of course heat pump’s specification. 
On the other hand, delta temperature occurred delivers additional 
thermal deformation and thermal stress that change the behaviour 
of material, especially friction mobilised at interface soil-concrete 
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piles also head uplift displacement and axial forces at interface 
concrete piles-upper structure. 

A preliminary numerical simulation is conducted with finite 
difference method under static thermal load for only cooling and 
only heating operation. The study considered two different condi- 
tions of contact between soil and an energy pile: perfectly contact 
and sliding contact with frictional interface elements. In this case, 
the pile is freely restrained both at its ends. The results show that 
by the presence of interfaces at zone of contact between soil and 
pile, the stresses and displacements obtained are lower than those 
with perfect contact. Cooling cycle in winter mode results in pile’s 
contraction so that produces surface’s settlement and temperature- 
induced tensile forces at pile. If interfaces are imposed at surface of 
contact, cooling cycle increases frictions mobilised at the upper part 
of pile while reduces at the lower part of pile. Reverses responses 
are obtained at heating cycle in summer mode. A higher gradient 
temperature applied will possibly create negative friction at inter- 
face soil-structure. There is a node at the middle part of pile that is 
not influenced by temperature-induced friction mobilised, which is 
called null point. This behaviour of energy pile under thermal load 
is none of difference with the results of in situ tests at Lausanne 
and Lambeth College. 

Numerical simulation under thermal static loads for only cool- 
ing/heating operation pointed out that there is slight influence in 
mechanical behaviour of pile. The problem in this kind of operation 
system lies in quantity of energy balance. Mechanical behaviour 
will be affected in seasonal operation time since the piles are sub- 
jected by cyclic thermal loadings during annual operation. It is 
suggested to concentrate on behaviour of interface due to cyclic 
solicitation in thermal regime with appropriate cyclic behaviour 
law of interface soil-structure considering coupling with temper- 
ature function for further studies. 
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